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摘 要 : 针对 传统 的 基于 稀疏 表示 的 DOA 估计 算法 单纯 利用 信号 
进行 稀 0 随 着 
结构 特征 和 统计 特征 , 提出 ee 下 DOA 


影响 信号 各 玻 重 构 效果 的 问题 。 使 用 分 块 大 疏 理论 对 信和 号 
往往 呈现 目标 群 测 向 的 特点 ， 为 了 能 够 更 好 地 利用 信号 的 


估计 算法 ， 使 用 块 稀 朴 理论 挖 据 信 号 的 内 部 结构 ， 充 分 利用 了 信和 号 


局 
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空域 稀疏 性 ， 寻 致 在 低 信 唆 比 时 稀 芍 性 能 变 差 ， 


的 块 内 稀疏 


目标 增多 及 作战 任务 改 1 DOA 估计 


性 和 块 间 相 关 性 ， 提 高 稀疏 重 构 性 能 ， 


进而 对 DOA 估计 效果 有 很 大 的 提升 。 仿 真实 验 表明 ， 相 比 于 经 典 的 DOA 方法 ， er 
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DOA estimation based on improved block sparse Bayesian learning algorithm 
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Abstract: The traditional DOA estimation algorithm based on sparse representation only uses sparse signal in spatial 


domain, which leads to poor sparse performance at low SNR and affects the effect of sparse signal reconstruction. This 


paper used block sparse theory to decompose the signal sparsely. With the increase of target and the change of combat 


mission, DOA estimation often presents the characteristics of target group direction finding. In order to make better use of 


the structural and statistical characteristics of signals, this paper proposed a block sparse DOA estimation algorithm based 


on Space-time combination (STC-BSBL). This method improves the performance of sparse reconstruction, and thus greatly 


improves the DOA estimation effect. Simulation results show that the proposed method has better estimation performance 


compared with the traditional DOA method. 
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前 ， 基 于 稀 疏 表示 的 雷达 信号 波 达 方向 (direction of 
arrival, DOA) 估 计 估 计算 法 广 受 关注 , 李子 高 等 人 利用 目前 
广泛 应 用 的 信号 稀 琉 表示 ， 使 用 是 优化 类 的 正 交 匹配 跟踪 
(orthogonal matching pursuit，OMP) 算 法 实现 信号 的 DOA 估 
计 ， 但 这 种 方法 只 考虑 了 信号 整体 的 稀疏 性 ， 往 往 会 造成 低 
信 噪 比 条 件 下 难以 精确 重 构 信 号 的 问题 。 实 际 上 ， 除 了 一 般 
稀 玻 性 ， 块 稀 玻 结构 在 DOA 估计 中 广 - 泛 存在 ， 如 今 的 测 向 
问题 ， 空 域 中 感 兴趣 的 目标 信号 往往 以 集群 的 形式 出 现 ， 群 
标 数量 众多 ， 和 群 内 分 布 密集 ， 呈 现 分 块 稀疏 性 。 文 献 [2] 
提出 块 稀 疏 分 解 的 信号 重 构 方法 ， 利 用 信和 号 块 的 结构 特性 提 
高 重 构 性 能 。 张 智 林 等 人 局 将 这 一 理论 进一步 发 展 ， 提 出 了 
块 稀 疏 贝 叶 斯 学 习 算 法 (block sparse Bayesian 
BSBL)， 创 新 性 地 利用 了 信号 的 分 块 稀疏 性 和 块 
现 了 精确 的 稀 朴 信号 重 构 。 
本 文 以 BSBL 为 基础 ， 将 其 拓展 到 多 测量 矢量 (MMV) 模 

型 DOA 估计 的 应 用 中 去 ， 在 考虑 信号 分 块 稀疏 性 和 信和 号 相 
关 性 的 基础 上 ， 兼 顾 多 观测 矢量 间 的 时 域 相关 信息 ， 以 窄带 


人 


油 


9 


learning ， 


间 相 关 性 实 


言 号 为 例 ， 提 出 使 用 空 时 联合 块 稀疏 贝 叶 斯 学 习 (sSpace-time 
combination BSBL, STC-BSBL) 方 法 对 其 进行 DOA 估计 。 

1 ”阵列 结构 及 数据 模型 

1.1 MMYV 模型 


M 个 阵 元 组 成 的 均匀 线 


考虑 工 个 远 场 窄带 信号 入 射 到 | 


收 稿 日 期 : 2018-08-15; 修 回 日 期 : 2018-10-10 
作者 简介 : 王 书 豪 〈1994-)， 男 ， 山 东 泰 安 人 ， 
肥 人 人， 教授， 博时， 主要 研究 方向 为 信号 与 信息 处 理 ， 


硕士 研究 生 ， 主 要 研究 方向 为 信号 与 信 ， 


d， 阵 列 接收 模型 


假设 信号 的 


间距 取 为 4= 和 2 ， 


阵 上 ， 一 元 线 阵 沿 轴线 均匀 放 


t， 每 个 阵 元 间距 相同 ， 均 为 


如 图 1 所 示 。 


1 阵列 结构 示意 图 


Fig. 1 Array structure 


P 心 频率 为 1 ， 带 宽 为 B， 波 长 为 和， 阵 元 
入 射 角度 为 90={9,8…9.} ， 则 t 时 刻 第 m 


个 阵 元 的 接收 数 : 


~ 3 


居 表 示 为 
ya (tl) = -Ti) +n, (0) 


m=1,2,...M (1) 


On Dasin(0),m=1,2,...M,i=1,.…,L 。 


如 果 考 虑 整个 天 线 阵 列 中 M 个 阵 元 接收 到 的 信号 , 将 接 


息 处 理 《540675575@qq.com); 
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收 数据 用 矩阵 的 形式 表示 为 ysvec(Y”) 
ry [eipnsin@) oipnsin(@).,. 0-jBusin(O.) x(D) nt) - A | 
0) _ Oiasin(0) -jlsntb) ,,, o-Psin(0,) 00) 六 人 四 xs 人 vec(X 7)J(vec(X 
: : : : 1 人 vec(N ) 
LDO Lee 其 中 xeRW ， AeR”*"，yeR"" 。 得 到 单 观 测 矢量 下 的 
将 上 述 过 程 拓展 到 多 观测 矢量 (multi measurement vector，” 块 稀疏 贝 叶 斯 学 习 的 DOA 估计 模型 
MMV) 模 型 ， 设 定 快 拍 数 为 P， 即 对 空域 内 所 有 感 兴趣 的 信 y=Axtn (7) 
号 进行 P 次 观测 ， 得 到 了 个 如 式 (2) 所 示 的 观测 矢量 ,用 和 矩阵 在 这 一 模型 下 ,信号 x 是 分 块 稀 疏 的 ,每 一 个 信号 块 二 大 
的 形式 表示 为 小 为 dxP， 其 信号 模型 示意 图 如 图 2 所 示 。 
Y=®DX+N. (3) FD 
YeR”* 是 P 个 测量 矢量 组 成 的 观测 矩阵 ，@ 是 由 空域 
密集 网 格 划 分 构成 的 过 完备 基 和 矩阵 ， 在 MMYV 模型 中 ， 称 为 [2 N] 
感知 矩阵 。 待 重 构 信 号 XeRY ， N, 为 加 性 高 斯 噪声 。 本 文 
进行 DOA 估计 的 关键 就 是 求解 信号 源 矩 阵 X 中 不 同 列 向 量 
之 间 非 零 信 号 的 位 置 ， 即 对 应 信号 的 方位 角 。 
1.2 块 稀疏 MMV 模型 
在 上 述 MMYV 模型 中 ， 本 文 作 以 下 假设 : 信号 x 中 不 同 
列 向 量 间 非 零 位 置 相同 。 通 过 这 一 假设 ，X 就 具备 了 行 稀疏 
(row sparse) 性 内， 非 零 行 矢量 具备 了 分 块 稀 玻 的 基础 ， 块 稀 图 2 ” 空 时 块 稀疏 信号 示意 图 
疏 的 优势 在 于 它 能 充分 利用 信号 内 部 的 结构 特性 ， 以 及 非 零 Fig.2 Space-time block sparse signal diagram 
块 之 间 的 时 域 相 关 性 ， 提 高 信号 的 稀疏 重 构 性 能 ， 应 用 到 正如 图 1 中 阵列 接收 ， 远 场 信号 往往 呈现 多 目标 、 集 群 
DOA 估计 中 ， 可 以 得 到 更 好 的 估计 效果 ， 更 快 的 收敛 速度 。 化 的 特点 ， 这 就 造成 了 目标 信号 不 仅 在 空域 (MMYV 模型 的 
下 面 建立 与 之 对 应 的 块 稀疏 MMYV 模型 。 将 信号 矩阵 XX  ” 行 向 量 ) 是 分 块 稀疏 的 ， 而 且 时 域 回 波 (MMYV 模型 的 列 向 


分 为 g 块 ， 非 零 信 号 块 个 数 W 满足 w<s， 第 P 个 信号 块 用 量 ) 也 存在 很 强 的 相关 性 中。 如 何 利用 这 两 类 信息 ， 研 究 快 
x 表示， 每 个 块 大 小 为 4;，d; 大 小 可 以 不 相等 ， 这 一 点 区 别 。 速 、 高 性 能 重 构 算法 , 是 本 文 最 关心 的 问题 。 基于 以 上 考虑 ， 


于 传统 块 稀疏 重 构 算 法 。 其 分 块 结构 式 如 下 本 文 提出 了 空 时 联合 块 稀疏 贝 叶 斯 学 习 算 法 ， 在 BSBL 算法 
X =[X, ,x NaH Kata sa Xy] (4) 的 基础 上 ， 将 信 号 的 分 块 稀 朴 JIL1 特性 ciock Sparse )、 块 内 相关 
性 结构 ee 0 以 及 多 观测 矢量 时 域 相关 性 


结构 联合 5 重 能 。 


,为 对 应 第 i 个 信号 块 的 Mxd 维 基 矢 量 。 引入 超 参量 xy, ，B ， 其 中 ，y 代 表 第 i 行 向 量 % 与 观测 
a SD 矩阵 Y 的 相关 性 , 信号 的 x 较 大 , 噪声 的 yx 较 小 , 在 参数 yx 的 

bd: 会 所 Im 学 二 ; ， 史 
2 空 时 联合 的 块 稀 玻 只 时 斯 学 习 算 法 学 习 中 ， 其 值 越 大 ，x 就 越 有 可 能 是 信号 分 量 ， 其 值 越 小 ， 
本 ee 块 稀疏 贝 叶 斯 学 习 框架 下 越 有 可 能 是 噪声 分 量 ， 一 旦 确定 了 x 中 的 非 零 位 置 ， 信 号 的 


的 信号 DOA 估计 模型 ， 然 后 给 出 相应 的 算法 步 又。 首先 建 方位 角 随 之 确定 。 因 此 对 超 参 量 的 学 习 是 核心 。 而 对 于 块 
立 贝 叶 斯 框架 下 块 稀 琉 模型 ， 然后 针对 多 测量 条 件 下 群 目标 稀 玻 信号 而 言 ， 还 需要 考虑 块 内 结构 ， 引 入 参量 6 ， 定 义 
信号 所 具有 的 空 时 相关 性 后 ， 提 出 空 时 联合 稀 玻 的 思想 ， 推 A=7XpB 为 % 的 协 方差 窍 阵 ，4 也 是 一 个 正定 对 称 矩 阵 ， 用 来 
导 了 模型 中 时 域 特 性 参数 和 结构 特性 参数 的 迭代 求解 方法 ， 描述 信号 X 的 空域 相关 性 。 
实现 了 信号 的 重 构 ， 进 而 实现 对 目标 信号 源 的 DOA 估计 。 对 于 时 域 相关 性 的 建 模 ， 本 文 从 实际 情况 出 发 ， 假 设 信 
2.1 贝 叶 斯 框架 下 块 稀疏 模型 的 构建 号 矩阵 XX 的 列 向 量 都 满足 同一 个 时 域 相关 结构 B ， 则 
对 块 稀疏 模型 的 处 理 常 使 用 稀疏 贝 叶 斯 学 习 的 方法 ， 这 p(%;A,B)= N(x;0,4 ®B) (9) 
是 因为 稀 玻 贝 叶 斯 学 习 不 但 可 以 方便 地 对 各 类 物理 信息 进行 进一步 假设 信号 矩阵 X 中 各 行 向 量 不 相关 , 则 式 (9) 可 以 
建 模 描述 ， 而 且 其 重 构 性 能 优异 且 对 相关 性 较 强 的 感知 矩阵 写 为 
有 较 好 的 鲁 棒 性 四。 稀 疏 贝 叶 斯 学 习 (sparse Bayesian p(XxIf,B)= N(xX;0,T ®B) (10) 
learning，SBIL) 由 Tipping 中 首次 提出 ， 并 经 过 Wipf 等 人 不 断 T=diag™ ({A, As,… An}) 
改进 ,在 稀疏 信号 重 构 领域 收 到 了 越 来 越 多 研究 学 者 的 青睐 。 假设 噪声 服从 零 均 值 高 斯 分 布 ， 即 
该 方法 通过 对 观测 y 和 信号 x 建立 分 层 贝 叶 斯 模型 ， 采 用 人参 No ~ N(0,4) vp,q 
数 化 的 概率 模型 型 对 系统 求 解 ， 为 求解 稀疏 的 结构 信息 提供 了 则 py DD ~ NC. aD daD 
种 非常 灵活 的 数学 方法 
p(y|B)~ N(y; PX, Bp"1) (5) 新 的 感知 矩阵 Ye@1/ 过 于 巨大 (其 大 小 为 PM xPN )， 
PSN)~ N(x;0,7) i=bN 在 稀 玻 重 构 算法 中 ， 存 储 和 计算 每 一 个 感知 矩阵 块 都 要 耗费 
于 这 一 方法 是 在 单 观测 矢量 (SMV) 的 前 提 下 提出 的 ， 很 大 的 存储 空间 ， 因 此 ， 考 虑 引入 矩形 正 态 分 布 ,将 SMV 


寻 此 ， 本 文 把 上 文 所 建立 的 块 稀 跑 MMYV 模型 转换 为 单 观测 。 ”模型 转换 为 矩阵 形式 下 的 贝 叶 斯 学 习 模型 101。 

矢量 下 的 块 稀 疏 贝 叶 斯 模型 ， 从 而 套用 上 式 进行 分 析 。 此 更 新 式 (10)(11)， 得 到 MMYV 模型 下 X 的 先 验 分 布 
利用 Kronecker 积 中 将 多 测量 矢量 映射 为 单 测量 矢量 ， 和 观测 信号 模型 

体 步 又 为 DCXIr,B)=MN(X:0.F.B) (12) 
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p(Y|4,B) = MN(Y: PX, 41,B) (13) 
结合 先 验 概率 模型 式 (12)(13), 得 到 似 然 函数 和 后 验 概率 
公式 为 
Pp(YI4,7,B)=p(Y|X,4) p(XIT,B) 
= MN(Y:0,C,B) (14) 
p(XIY;4,5,B)=N(p,>) (15) 
其 中 
C=A/1+ PTO 
均值 为 
p= /3p'Y 


2=(T'+Y p09) 
首先 ， 对 超 参 量 4 的 估计 ， 常 把 它 当 做 一 个 罚 参 数 ， 根 
据 信 品 比 的 不 同 ， 对 4 赋予 不 同 的 值 5 


0.1lyl, SNR < 20dB 
4=40.01|yl, SNR < 20dB 
10™, 无 噪声 


对 均值 和 协 方差 的 求解 等 价 于 对 超 参量 4、8 的 估计 问 
题 ， 正 确 估计 出 超 参 量 即 可 获得 最 大 后 验 概率 的 解 ， 从 而 重 
构 X， 得 到 所 需 的 DOA 估计 ， 为 此 ， 构 造 如 下 的 对 数 似 然 
函数 : 


L=M log|lBl+ Plog|C|+7r[B"Y'C''Y | (16) 


通过 对 似 然 函 数 中 各 个 超 参量 求 偏 导 ， 可 得 到 各 参量 的 
更 新 公式 , 为 了 简化 计算 , 将 上 式 写成 两 个 代价 函数 的 形式 。 


L(B)s M1og|lBI+Tr[B"'Y'C''Y | (17) 


L(A)sPlog|Cl+7Tr[B"'Y"C"Y] (18) 


在 ZL(B) 中， 将 驴 看 做 已 知 变 量 ，ZL(4) 中 ， 将 8 看 做 已 
知 变量 , 通过 给 8 一 个 初始 化 的 值 ,采用 式 (17) 与 式 (18) 对 超 
参量 进行 交 蔡 更 新 M1， 步骤 为 ， 首 先 初 始 化 B=/ ， 更 新 (4) 
得 到 A， 随后 计算 得 到 C， 继 续 对 ZL(B) 更 新 得 到 新 的 B。 经 
过 次 迭代 后 ， 驴 ，8B 收敛 于 一 稳定 值 。 

对 于 参量 8 的 更 新 ， 可 以 直接 对 ZL(B) 求 偏 导 


oL(B) _ 1_ pivriAl 1 
TE MB-B'YC'YB (19) 
令 22 到 -0， 得 
YI7CTY 
B= (20) 
对 块 相关 算 阵 4 的 更 新 ， 代 价 函 数 5(4) 中 的 C 可 写 为 
C=D +bA97 (21) 


ee A T 
| 中 定义 D,=AI+ P09, 


mi 


应 用 和 矩阵 求 逆 定理 [器 并 对 代价 函数 进行 最 小 化 处 理 可 


en 


ss9D'Y 
qs D0 


,B'g” 
A -js C2) 
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通过 计算 得 到 的 4 ， 通 过 x= 寺 六 [A] 计算 得 到 上 9。 估 
计 出 所 有 的 超 参量 后 ， 后 验 概率 密度 得 出 ， 即 


X=p=(42" +D PPY 


= (NU+9 39P) Y (23) 

从 而 恢复 出 原 信 号 ， 并 根据 x 的 非 零 位 置 ， 完 成 信号 的 

DOA 估计 。 使 用 流程 图 来 描述 整个 信号 重 构 及 DOA 估计 过 
程 ， 如 图 3 所 示 。 


初始 化 4 ，A，，B 的 值 


i 


一 一 > ”计算 均值 k 和 方差 二 


7%，B 收 敛 于 一 稳定 值 ? 


计算 后 验 概率 X* = 


根据 XX* 恢 复原 信和 号， 
完成 得 到 DOA 佑 计 
图 3 基于 STC-BSBL 算法 的 DOA 估计 流程 图 
Fig.3 DOA estimation flow chart based on STC-BSBL algorithm 


3 ”仿真 实验 


3.1 DOA 估计 效果 比较 

设 定 均匀 线 阵 阵 元 数 为 30, 阵 元 间距 为 信号 频率 对 应 的 
半 波 长 ， 空 域 范围 [-90",90-] ， 真 实 信号 数目 3， 真实 信号 来 波 
方位 角 9=[0,35 ,40]i=123 ， 输 入 信 品 比 SNR=8dB， 快 拍 数 
P=100。 进 行 1000 次 蒙特 卡 洛 实验 ，DOA 估计 结果 如 图 4 
所 示 。 

将 凸 优 化 类 正 交 匹配 追踪 (OMP) 算法 和 基于 子 空间 类 
的 MUSIC 算法 对 信号 进行 DOA 估计 ， 与 本 文 所 使 用 的 
STC-BSBL 算法 的 DOA 估计 效果 进行 比较 。 信 和 号 数目 为 4， 
来 波 方向 分 别 为 6=[30 ,45 ,60.75] ， 其 他 条 件 同 上 文 不 变 。 从 
仿真 结果 如 图 5 所 示 。 


2.5 


T T T T 
一 + 一 Estimation of DOA 
口 The real DOA 


wb 
n 


Amplitude(dB) 


Angle(®) 


妈 4 ”STC-BSBL 算法 的 DOA 估计 效果 
Fig.4 DOA estimation effect of STC-BSBL algorithm 
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图 5 不 同 算法 的 DOA 估计 效果 对 比 

Fig.5 Effect Comparison of different algorithms of DOA estimation 
从 仿真 结果 可 以 看 出 ， 在 给 定 的 条 件 下 ， 三 种 方法 均 能 

实现 对 信号 的 DOA 估计 ， 但 从 估计 精度 的 角度 上 看 ， 

STC-BSBL 算法 和 OMP 算法 拥有 比 MUSIC 算法 更 罕 的 角度 

估计 范围 ， 和 更 低 的 幅 值 适 应 度 。 

3.2 DOA 估计 正确 率 对 比 

目前 ， 基 于 稀疏 分 解 的 DOA 估计 算法 除了 凸 优化 类 算 

法 ， 还 包括 稀疏 贝 叶 斯 学 习 框架 下 常用 的 块 稀 朴 贝 叶 斯 学 习 

算法 (BSBL)， 基 于 时 序 结构 的 稀疏 贝 叶 斯 学 习 (TM-SBL) 

05, 本 实验 以 稀 玻 度 K 为 自 变量 , 仿真 了 这 三 种 方法 的 DOA 

估计 正确 率 并 与 本 文 方法 进行 对 比 。 四 种 方法 在 不 同 稀 玻 度 

下 的 识别 正确 率 曲线 如 图 6 所 示 。 


人 放生 一 看 


STC-BSBL 
OMP 中 
一 ~ 一 TM-SBL 
一 ”一 BSBL 中 


09 上 + 


0.8r 


0.7 上 


0.6 上 


0.5 


估计 正确 率 


0.4 上 


0 5 本 本 15 20 25 
图 6 识别 正确 率 
Fig.6 The correct recognition rate 

从 仿真 结果 可 以 看 出 ， 本 文 使 用 的 STC-BSBL 算法 在 同 
等 的 测量 条 件 下 有 着 更 好 的 重 构 性 能 ， 在 同等 稀疏 度 下 ， 有 
着 更 高 的 DOA 估计 正确 率 ， 且 能 适应 更 多 数目 的 信号 DOA 
估计 。 
3.3 仿真 实验 3 DOA 估计 均 方 根 误差 分 析 

为 了 能 更 好 地 凸显 本 文 算法 在 DOA 估计 的 优越 性 ， 分 
别 以 信 品 比 和 快 拍 数 为 自 变 量 ， 比 较 五 种 方法 在 进行 DOA 
估计 是 的 均 方 误差 曲线 。 定 义 DOA 估计 的 均 方 根 误 差 为 


RMSE=| 二 > > (6, -4 (24) 


结果 如 图 7、8 所 示 ， 从 图 中 可 以 看 出 ， 随 着 信 噪 比 和 快 
拍 数 的 提高 ， 每 种 算法 的 RMSE 均 在 减 小 ， 但 在 低 信 噪 比 和 
小 快 拍 数 时 ，OMP 算法 和 ISSM 算法 的 估计 误差 明显 增 大 。 
BSBL 和 TM-SBL 算法 利用 了 信号 的 块 内 结构 或 时 域 相 关 结 
构 ， 重 构 误 差 有 一 定 降低 ,但 STC-BSBL 算法 由 于 联合 了 信 
号 的 空域 相关 性 和 时 域 相关 性 ， 估 计 误 差 始 终 是 最 小 的 。 
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RMSE() 


SNR(dB) 
图 7 均 方 根 误差 随 信 噪 比 的 变化 曲线 
Fig.7 ， Variation of RMSE with the SNR 


STC-BSBL 
一 OMP 
一 一 TM-SBL 
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一 扩 ISSM 


RMSEO) 


快 拍 数 
图 8 均 方 根 误差 随 快 拍 数 的 变化 曲线 


Fig.8 Variation of RMSE with the number of snapshot 
4 结束语 


本 文 在 稀疏 贝 叶 斯 学 习 框架 下 ， 分 析 了 块 稀 疏 理论 在 稀 
朴 信 号 重 构 方 面 的 优越 性 ， 推 导 了 块 稀疏 贝 叶 斯 学 习 模 型 ， 
并 针对 多 观测 矢量 模型 下 的 窄带 信号 DOA 估计 问题 ， 构 建 
了 块 稀 疏 MMYV 模型 ， 在 现 有 的 BSBL 算法 的 基础 上 进行 改 
进 ， 在 考虑 信号 块 内 相关 性 的 基础 上 ， 联 合 了 多 次 观测 下 信 
号 的 时 域 相 关 性 ， 提 出 空 时 联合 块 稀疏 贝 叶 斯 学 习 算法 。 通 
过 比较 实验 ， 证 明了 该 算法 对 雷达 信号 DOA 估计 在 测 向 精 
度 上 有 着 更 好 的 效果 。 
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